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Abstract 

'> . 

qq , Basic concepts of quantum integrable systems (QIS) are presented stressing on the 



unifying structures underlying such diverse models. Variety of ultralocal and nonul- 
tralocal models is shown to be described by a few basic relations defining novel algebraic 
entries. Such properties can generate and classify integrable models systematically and 
also help to solve exactly their eigenvalue problem in an almost model-independent way. 
The unifying thread stretches also beyond the QIS to establish its deep connections with 
statistical models, conformal field theory etc. as well as with abstract mathematical 
objects like quantum group, braided or quadratic algebra. 



1 Introduction 

A number of quantum models of diverse nature shows the important property of integra- 
bility in low dimensions, enabling us to solve them exactly. Examples of such models may 
extend from field models like sine-Gordon (SG), nonlinear Schrodinger equation (NLS), 
derivative NLS (DNLS) etc. to the discrete models like Toda chain (TC), relativistic TC, 
isotropic (XXX) or anisotropic spin chains (XXZ,XYZ) etc. Similarly they may have 
varied basic commutation relations ranging from ultralocal to nonultralocal models. How- 
ever, the fascinating feature of such models, as revealed in recent years, is their common 
underlying algebraic structures defined through a basic relation called Yang-Baxter equa- 
tion. The unifying spirit goes also beyond the domain of the quantum integrable systems 
(QIS) to reveal its deep relations with variety of other disciplines, which were thought to 
be completely unrelated in the recent past. Thus in one hand the QIS is deeply connected 
with abstract mathematical objects like quantum group, noncocommutative Hopf algebras, 
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quadratic and braided algebra, universal 7£-matrix etc. and on the other hand it has in- 
timate relation with stat-mech problems, conformal field theory (CFT), knot and braid 
theories etc. 

The starting of QlS-age should possibly be counted from 1931, when Hans Bethe was 
the first to solve exactly an interacting many-body quantum problem, e.g. the isotropic 
Heisenberg spin chain, through his by now celebrated ansatz called the Bethe ansatz fij. 
This was followed by many other applications and successes of the method in solving var- 
ious physically interesting models || like XX Z-spin chain H|, many-body problem with 
attractive j|] or repulsive || 5-function interactions, Hubbard model §, Kondo problem 
Q etc. Next came the landmark discovery of Rodney Baxter, in which he not only solved 
the completely anisotropic XYZ spin chain but also established its deep link with a stat- 
mech problem, e.g. the eight-vertex model Q. However, the genuine theory of QIS based 
on all the preceding achievements in quantum as well as classical integrable theories was 
started taking shape only in the late seventies || fl(| . Synthesising the concept of the Lax- 
operator, borrowed from the classical theory, along with the Bethe ansatz the quantum 
inverse scattering method (QISM) was developed mainly through the works of the Russian 
school, where the celebrated Yang-Baxter equation plays the central role and an algebraic 
formulation of the Bethe ansatz was invented for covering wider range of models including 
field theoretic models [11]. The connection of the integrable system with the conformal 
field theory (CFT) |l^] [ 13 1 as well as with the quantum group [|l4|] and the braid and knot 
theory |15| was understood quite recently. Thanks to this ever-growing progress, today the 
theory of QIS is considered to be a major and important branch of theoretical and math- 
ematical physics, with many research groups all over the globe engaged in active research 
in the field. 

We consider the notion of integrability in the Liouville sense, i.e. call a system integrable, 
when the number of independent conserved quantities coincides with the degree of freedom 
of the system. In the quantum case they correspond to conserved operators with one being 
the Hamiltonian, commuting with them. Moreover, their independence demands that all of 
them must form a mutually commuting set of operators, which in field models are infinite 
in number. We will see that the quantum Yang-Baxter equation, which plays the role 
of the basic equation in QIS, ensures in fact such a property. Though in physics we are 
interested mainly in the Hamiltonian and its eigenvalue solutions, a completely integrable 
system actually guarantees much more, it should in principal give the whole hierarchy of 
conserved operators along with their exact eigenvalue solutions, which is one of the main 
aims of the QISM. 

We present here the basic concepts of QISM with a major stress on the unifying scheme 
for both ultralocal and nonultralocal quantum models. We elaborate on the fascinating and 
novel algebraic structures revealed by the theory, which help to generate and solve such 
models in a systematic way. We also focus on the intriguing connections of the theory with 
the stat-mech and CFT models. The section headings explain the topics dealt in each of 
them. 
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2 Lax operator and examples of integrable systems 



The central idea of classical ISM [16] is that, instead of dealing with the nonlinear equation 



in (1 + l)-dimensions directly, it constructs the corresponding linear scattering problem 

%{x,X) = L(u{x,t),p(x,t),X) T(x,A), (2.1) 

where the Lax operator L(u,p, A) depending on the fields u,p and the spectral parameter 
A contains all information about the original nonlinear system and may serve therefore as 
the representative of a concrete model. The field u acts as the scattering potential. The 
aim of ISM is to solve the inverse problem by first finding a canonical mapping from the 
spectral data, which are like action-angle variables, to the scattering potential, i.e. the 
original field and using it to construct the exact solution for the given nonlinear equation. 
Soliton is a special solution, which corresponds to the reflection-less potential. 

The QISM succeeds to generalise the notion of the Lax operator also to the quantum 
case. Since the canonical variables u,p or etc. now become operators acting on 

the Hilbert space, the quantum L(A)-operators are unusual matrices with non-commuting 
matrix elements. This intriguing feature leads to nontrivial underlying algebraic structures 
in QIS. We present below a few concrete examples of the Lax operators associated with well 
known models to give an idea about the structure of this immensely important object in the 
integrable theory. The field models are represented by field equations and the corresponding 
continuum Lax operator C, while the lattice models are given through the Hamiltonian H 
and the discrete Lax operator L. 

I. Trigonometric Class: 

1. Sine-Gordon (SG) model 



/ \ f \ m2 • t i w /- I Wi msin(A — rru) \ 

u(x, t)u - u(x, t) xx = sm(f]u(x,t)), Lsg = . N . , p = u 

rj \ msm{A + r]U), —ip I 



2. Liouville model (LM) 



3. Anisotropic XXZ spin chain 



N 

H = XX^n+l + ^n+l + cos wlol+i ) , 



sin(A+ §cr£), sinr?o- n 
sin 7/ct+, sin(A-2cr^, 



II. Rational Class: 



(2.2) 



L(x,t) tt - u(x,t) xx = I e 2 ^), C LM = i ( j • (2.3) 



L n (0 = i I 2 nh . ,/ (2.4) 
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1. Isotropic XXX spin chain 

N 3 t A + io* 



H = J2 E(« +1 ), L n (0 =i( i~ 2 \ 3" (2.5) 
2. Nonlinear Schrodinger equation (NLS) 



%il>(x,t) t + ^(x,t) xx + V (^(x.t)iP{x,t))^{x,t) =0, £tvl5(A)=( ?' ^ ]. (2.6) 
3. Toda chain (TC) 

H = J2 Qrf + . ^n(A) = ( ^ e \ ) . (2.7) 

Notice that we have put the models under different classes and apparently diverse looking 
models into the same class. The meaning of this will be understood as we go further, though 
it might have already sent a signal of the fascinating unifying feature in the integrable 
systems. Observing carefully we also see that the off-diagonal elements (as tp, ifi* in (p.6| ) 
and a~,cr + in (|2.4|)) involve creation and annihilation operators while the diagonal terms 
are the number like operators. It is crucial to note that under matrix multiplication this 
property is preserved, which as we will see below, has important consequences in their 
algebraic Bethe ansatz solution. 

3 Yang-Baxter equation, i?-matrix and notion of quantum 
integrability 

As we see, the Lax operators are local functions of x, or if we discretise the space of every 
lattice site i. However, since the integrability is defined through the conserved quantities, 
which are global objects, in proving integrability we evidently need some global entries 
constructed from the local Lax operators. Such an object can be formed by matrix multi- 
plication of the Lax operators at all lattice sites as 

N 




r(A)=nw= . (3.1 




As we have indicated above, the off-diagonal global operators B(X),C(\) are like cre- 
ation/annihilation operators and related to the angle- variables, while A(X), D(X) corre- 
spond to action variables. For ensuring integrability one must show that t(A) = trT(X) = 
A(X) + -D(A) generates the conserved operators : lnr(A) = J2jCjX J with [H,C m ] = 
and mutual commutation: [C n , C m ] = 0. This is in fact achieved by a key condition on the 
quantum Lax operators (for ultralocal models) given by 

R 12 (X - fi) L U (X) L 2l (p) = L 2i (p) L U {X) R 12 (X-fi). (3.2) 
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with the notations Lu = Li <S> I, L 2 { = I <8> Li. This matrix relation is known as the 
Quantum Yang-Baxter equation (QYBE), where apart from L-operator i?i2(A, ^)-matrix 
with c-number functions of spectral parameters and acting nontrivially only in the first two 
spaces V (8> V I appears, which in turn satisfies the YBE 

i?i 2 (A-/i) i?i 3 (A- 7 ) Rw(ji--y) = Rosin-i) Ri3(X--y) i*i 2 (A-/i). (3.3) 

Due to an additional ultralocal property of the L-operators: [Lu(X), L 2 j(fJ,)] = 0, i ^ j of a 
large class of models, one can treat these 'specially separated' operators almost as classical 
objects. Therefore multiplying them for all sites and repeatedly using ( |3.2D one may arrive 
at the same QYBE 

#12 (A - n) 2i(A) T 2 (ji) = T 2 (jj) Ti(\) R 12 (X - fi), (3.4) 



but for the global object T(A) defined in (|3.l|) . This invariance relation for the tensorial 
product reflects a deep algebraic property related to the Hopf algebra, to which we will 
return later. Taking now the trace of relation fl3.4|) , ( since under the trace i?-matrices can 
rotate cyclically and thus cancel out) one gets [r(A), r(/i)] = 0, proving the commutativity of 
C n for different re's. Thus starting from the local QYBE (3.2) and following several logical 



R(X) 



(3.5) 



steps we finally establish the complete integrability of the quantum system. Therefore 
the validity of (|3.2j ) may be considered to be sufficient for the quantum integrability of the 
ultralocal systems and the local QYBE as the basic equation in the QIS. We will show below 
that even for nonultralocal models one can formulate the local QYBE in some generalised 
form and prove the integrability, though the procedure is much more involved. The 4x4 
i?(A)-matrix solution of YBE |3.3| is rather easy to find, which in the simplest form may be 
given as 

( fW \ 

1 /i(A) 
A(A) 1 

V J 

There are usually only two types of solutions (we shall not speak here of more general 
elliptic solutions), namely trigonometric with 

/(A) = ^±^, /x(A) = ^f (3.6) 
sin A sin A 

and the rational with 

/(A) = ^, hW = J- (3-7) 
In the examples of integrable systems the trigonometric and rational classes mentioned are 
associated with the respective i?-matrices presented here. 



4 Universality in underlying algebraic structure 

In sect. 2 we have given the Lax operators for a number of models and grouped some like 
SG, LM and XXZ in the same trigonometric class, while equally varied models like NLS, 
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TC, XXX etc. are in the same rational class. Here we address to the intriguing question 
that why a wide range of diverse models share the same R{\) and construct also their Lax 
operators exploring the universality feature underlying such systems. 

To generate models of the trigonometric class we start with the 72j r j g -matrix solution 
( |3,6j ) and look for the the Lax operator of the corresponding ancestor model in a generalised 
form 

L t (X)=i( sin ( A + 7 ? s3 )( c ^ + C D -icos(A + r?s 3 )(c^ - c[), (2sin? ? )5" \ 
\ (2sinr/)5' + , sin(A — r]s 3 )(c2 + C2 ) — i cos(A + r/s 3 )(c^ — C2 ) /' 

The abstract operators s^,S^ describe an algebra given by the following relations 



[ S 3 ,5±] = ±S ± ,[S + ,S~] = -\ (M+[2 S 3 ] (? + i^-cos(27 ?S 3 ) ) . (4.2) 

£ \ sin tj j 

where \x} = q ~ q _i = SlSiB q, = e "? anci jy± are related to the central elements of 

L 1L i q—q 1 sin 77 ' ^ 

the algebra : ,i = 1, 2 as = cfc^ ± c^cj. Note that (^2) is a new type of algebra 
underlying such QIS. This is in fact a quadratic algebra and unlike the known Lie algebra 
or its deformations there are arbitrary multiplicative central elements in the RHS of 
(|4.2j), which can take trivial values. For different values of these elements one gets different 
kinds of the algebra along with different types of the Lax operators. Therefore starting from 
this universal structure of the ancestor model ( |4.1| ) and considering various representation 
of (4.2) with concrete choices for , one can generate systematically the Lax operators of 
diverse models belonging to this class. 

Remarkably, for the choice M~ = 0, M + = — 2 one gets the well known quantum algebra 
(QA) U q (su(2)) with the defining relations 

[s 3 ,S ± ] = ±S ± , [S + ,S-] = [2s\ (4.3) 

Let us consider first the models related directly to this QA with the Lax operators obtained 
from ( |4.1| ) for the choice = = — c\ = —c^ = 1. The fundamental representation 

S ± = ±a ± ,s 3 = ±a 3 (4.4) 

constructs clearly the Lax operator of the spin-i XXZ-chain (|2.4[). Another nontrivial 



example obtained through the bosonic representation ( with [linjPm] — V 

) of QA: 



A 

3 n = u n , S-=g(u n )e iA ^,S+ = (S-)i, (4.5) 



where g(u n ) = [1 + |m 2 A 2 cos 2r](u n + i)]2 yields an exact lattice version of the sine- 
Gordon model. The Lax operator of the model can be derived readily from (4.1) using 
(fL5p- At A — » one gets the SG field model with the Lax operator obtained as a x L n = 
I + ACsg(x) + O(A) reproducing (gg). 

We now go beyond QA and consider more general structure ( |4.2| ). Note that the choice 
c~l = 02 = A, Ci = 02 = giving M + = M~ = A 2 corresponds to an intriguing algebra 

-2s 3 

with [S' + ,S' _ 1 = iA—. — . A bosonic realisation of this algebra as 

L ' J sin?7 



S 



3 



'11 > 



S+ = /K)e^,S+ = (S-)t, (4.6) 
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with f(u) = ( 1 + A 2 e ir, ( 2u+in ^) 2 constructs again from the same (4.1) the Lax operator of 



an exact lattice Liouville model and at A — » the Liouville field model (|2.3| ). We have 
witnessed here an important consequence of the the underlying algebraic structure, from 
which through simple realisations we could generate exact lattice versions of SG and LM, 



which were originally found in an intuitive and much involved way [38] [37]. 

We have shown thus that the Lax operators of the models, at least for the trigonometric 
class presented in sect. 2, can be constructed systematically and the set of diverse models 
put under the same class are in fact descendants of the same ancestor model ( t4.1| ). We 
continue this procedure to generate other members of this class like the quantum Derivative 
NLS (DNLS), Ablowitz-Ladik model, relativistic Toda chain etc. as different realisations of 
the unifying algebra (4.2). For example, at cf = cf = 1, c{ = i-^q, c^ = ij^ corresponding 



to M + = — sin??, M = i-f cos 7/, the quadratic algebra (|4.2f ) after the substitution 

S + = kA, S~ = kA\ s 3 = N, k = -z(cotr?)3 
leads directly to the well known (/-oscillator algebra 

[A, At] = — cos(7 ? (2A^ + 1)), [A, N] = A, A\N] = -A\ (4.7) 

COST] 

and therefore to a novel integrable (/-oscillator model with explicit L-operator and the Rtrig 
-matrix ( |3.6| ). A realisation of the (/-oscillators in bosonic operators: [Vvi, = in 
the form 

An = (f )* ^nM^, N n = ^ n (4.8) 

derives an interesting exact lattice version of the quantum DNLS model, which at the 
liniA-^o yields the Lax operator 



£dnls 



(£) = i (-^£ 2 <7 3 + ^{x)ip{x)K + £(V> t (x)<7 + + ip(x)a , k = diag(K-,-K+) (4. 



of the DNLS field model, which is exactly solvable through the Bethe ansatz [2C]. A fusion 
of such DNLS model in turn can generate the well known massive Thirring model [35]. 

A different set of choice for the elements cf likewise results from ([h^) other algebras 
and related integrable models. The simplest choice = 0, cf = 1 leads to a light-cone SG 
model, while = 0, cj~ = — cf = 1 giving = with realisation 

s 3 = p s~ = —e u , S~ = —e u (4.10) 
sin 7] sin rj 



gives from the ancestor model (4.2) the Lax operator of a quantum relativistic Toda chain 



34]. It should also be noted that under a twisting transformation: 

KigW = F(9)R trig (X)F(9), with F(9) = e^ 3 ® 1 ' 1 ®^ (4.11) 

one gets a twisted trigonometric i?-matrix and a corresponding ^-parameter extension of 
the algebra (|4.2j). Proceeding as above and choosing properly 9 one gets yet another set of 
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integrable models, a prominent example of which is the Ablowitz-Ladik model [33] [35]. This 
demonstrates clearly the unifying feature of the algebraic structure in generating a wide 
range of descendant models from the Lax operator fl4.1[ ) all sharing the same trigonometric 
ii-matrix inherited from the ancestor model. 

Let us consider now the q — ► 1 limit when Rtrig —> Rrat giving the rational i? rai -matrix 
(pj), which can be expressed also as R ra t(^) = A + HP. The ancestor model reduces to the 
corresponding rational form 



L(XU= T \„ , " • (4-12) 




As a consequence of the QYBE ( p.2| ) the K operators satisfy the algebra given through the 
relations 

[#12, K 21 ] = {clKl - c?K 2 °), [Kl K° 2 ] = 

[Kf,K 12 ] = eiK 12 cl [Kl K 21 ] = -e^c?, (4.13) 

with ei = 1,62 = —1 and clc 2 are the central elements. Note that though it is much 
similar to a Lie algebra, it is in fact again a quadratic algebra with arbitrary multiplicative 
central elements. Since L(X) operator given in the general form ( [4.12 ) satisfies QYBE 



and associated with the rational i? m t-matrix, it may be taken as the ancestor model for 
generating integrable models belonging to the rational class. 

i 



A nontrivial choice for the central elements c? = Cn = 1 and the notation 



, Kl = -Kl = s 3 , K 12 = s+, K 21 = s- (4.14) 

recovers the su(2) spin algebra 

[s 3 , s ± ] = ±s ± , [s + , s~] = 2s 3 , (4.15) 

which is a Lie algebra. The integrable models like the XXX spin chain and the NLS model, 
our acquaintances from sect. 2, can be generated from the case (|4.14|) . In fact the simplest 
spin-^ representation through Pauli matrices s a = \a a reduces ( 4. 12| ) to the Lax operator 



of the XXX chain (2.5). On the other hand bosonic representation of the spin operators 



given by the Holstein-Primakov transformation 

s 3 = s - AV^s" = A^(2s - A^VOM, s + = (s _ ) f (4.16) 
leads (4.12) to the quantum integrable Lattice NLS model. The continuum limit takes it 



to the NLS field model Q. 

The general algebra fl4.13|) permitting trivial values for a® however allows more freedom 
for constructing models. Choosing c 2 = 0, = i for example one can set Kl == 0, 
reducing the algebra to 

[Kl K 12 ] = iK 12 , [Kl K 21 ] = -iK 21 [K 12 ,K 21 ] = . (4.17) 



S 



A bosonic representation of the generators: 



K»=p,K 12 = e- u , K 21 =e u , (4.18) 

as is evident from ( |4.12| ) yields the Lax operator of the Toda chain (^[7]). Another simple 
lattice NLS model can be obtained taking again the values as c 2 = 0, = —An 2 with the 
representation 

Kl = nA 2 cpip + 1, K$ = 1, K 12 = iA^/nip, K 21 = -iA^cp (4.19) 

where the operators ip, <fi are canonical. At the continuum limit the same standard NLS 
field model (^) is recovered. 

Finally we can consider the twisted rational i?f at (A)-matrix through transformation 
( [1.11 ) with the same twisting operator F(6). As a result one gets a ^-deformation of the 



algebra ( 4.13 ) in which cj's no longer remain central elements. We can construct again in a 
similar way the integrable models belonging to this class. Thus one can obtain the lattice 6- 
deformed NLS, ^-deformed Toda chain, Tamm-Dancoff (/-bosonic model or Dzhyaloshinsky- 



Moriya spin chain, which is nothing other than the ^-deformed XXX chain |3q| |31] 

Thus the algebraic structure in the QIS clearly plays a crucial role in unifying diverse 
models of the same class as descendants from the same ancestor model and at the same 
time realisations like (|4.5| ) gives a criterion for defining integrable nonlinearity as different 
nonlinear realisations of the underlying quantum algebras ( |4.lD or ( 4.13j ). We see below 



that this universality feature is also present curiously in their Bethe ansatz solution, which 
facilitates their exact eigenvalue solution in an almost model-independent way. 

5 Universality in Bethe ansatz solutions 

Coordinate formulation of the Bethe ansatz, though more effective for finding the energy 
spectrum :H \ m >= E m \ m > of concrete models, dependents heavily on the structure of 
the Hamiltonian and consequently lacks the uniform approach of its algebraic formulation. 
We therefore focus briefly only on the algebraic Bethe ansatz method for solving the general 
eigenvalue problem 

t(A) I m >= A m (A) | m > (5.1) 

and highlight its universal feature. 

Apart from the integrability condition, the QYBE (3.4) represents also a set of commu- 
tation relations between action and angle variables given in the matrix form. They can be 
obtained by inserting matrix T in the form fl3.1| ) and the i?(A)-matrix solution as ( |3.5| ) in 
the equation fl3.4|) . Such generalised commutation relations dictated by the QYBE are of 
the form 

A(X)B(fi) = f^-X)B( f j,)A(X)-f 1 ( f j,-X)B(X)A( fJ ,), (5.2) 
D{X)B(n) = f(X-n)B^)D(X)-f 1 ( l i-X)B(X)D( t i), (5.3) 
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together with the trivial commutations for [-A(A), A([j,)] = [B(X), B(fi)] = [D(X), D(fj,)] = 
[A(X), D(fjL)] = etc. It is now important to note that the o/ /-diagonal element B(X) acts 
like an creation operator (induced by the local creation operators of L(X) as argued above). 
Therefore the m-particle state | m > may be considered to be created by -B(Aj) acting m 
times on the pseudovacuum | >: 

| m >= B(X l )B(X 2 ) ■ ■ ■ B(X m ) | > (5.4) 



If one can solve the general problem (5.1), the eigenvalue problem for all 



1 d 

C n = -t-^t-tat(A) |a=o (5-5) 
n! oX 

can be obtained simultaneously by simply expanding A(A) as 

d\m>= A^ l (0)A~ 1 (0) \m>, C 2 \ m >= (A^(0)A~ 1 (0)) / | m > (5.6) 

etc. , where say the Hamiltonian H = C\. 

Evidently for solving ( [iTl] ) through the Bethe ansatz we have to drag r(A) = A(X)+D(X) 
through the string of Bx^s without spoiling their structures (and thereby preserving the 
eigenvector) and hit finally the pseudovacuum giving A(X) | >= a(A) | > and D(X) | 
>= /3(A) | >. 

Notice that for this purpose (5.2, 573|) coming from the QYBE are indeed the right kind 



of commutation relations but for the second terms in both the RHS, where the argument 
of B has changed: B — > B spoiling the structure of the eigenvector. However, if we put the 
sum of all such unwanted terms= 0, we should be able to achieve our goal. In field models 
such unwanted terms are however absent, while in lattice models they may be removed by 
the Bethe equations induced by the periodic boundary condition giving 

This in turn serves as the determining equations for parameters Xj. Ignoring therefore the 
second terms for the time being and making use of the first terms only, as argued above we 
finally solve the eigenvalue problem to yield 

m m 

A m (a) = n f( x i - a )«( a ) + n /( a - ( 5 - § ) 

3=1 3=1 

The structure of the eigenvalue A m (A) reveals the curious fact about the Bethe ansatz 
result, that apart from the a(a),/3(a) factors the eigenvalue (|5.8| ) and the Bethe equation 
depend mainly on the nature of the function /(A — Xj), which are universal within 



a class of models and determined only by the choice of i?-matrix as (3.6) or (3.7). The 
coefficients a(A),/3(A) determined by the concrete form of the Lax operators and the def- 
inition of the pseudovacuum are the only model-dependent part. Therefore models like 
the DNLS, SG, Liouville and the XXZ chain belonging to the trigonometric class share 
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similar type of eigenvalue relations with individual differences only in the form of a(A) 
and /3(A) coefficients. Thus this deep rooted universality in integrable systems helps to 
solve the eigenvalue problem for the whole class of models and for the full hierarchy of 
their conserved currents in a systematic way. For example, for the XXZ model using the 
Lax operator structure (|2.4|) and | > defined as the all spins up state, one easily finds 



a(A) = sin 7V (A + 77), /3(A) = sin^ A. Using subsequently the i? tr j 9 -matrix information ([T(]) 
one derives from ( |5.7|) (with a shift A — > A + 2 ) the Bethe equation 

/ sin(A fc + | ) \ N _ " sin(A fc - A,- + rj) 

\ s ^{\ k -D) AJL s in(A fe - \j - rj) ' { ' } 
for j = 1, 2, . . . , m. Similarly ( p\q ) gives the eigenvalue 

N„ ^ sin(A 7 -A + f) "L s in(A - A,- + 32) / 

yielding for H xxz = Ci, the energy spectrum 



rn 



E<g}=A(\)'A- 1 (\) | A=0 = s in77^— +ATcotr ? . (5.11) 

i=1 sm(Aj - 2 J sm^Aj + 2 J 

At the limit r/ — > 0, sin A — * A, when the i?-matrix along with its associated models 
reduce to the rational class, one can derive the Bethe ansatz result of the isotropic XXX 
chain directly by taking the limits of the above expressions for XXZ model. The Bethe 
ansatz result for the NLS model of the same rational class also have very similar structure. 

It should be remarked that unlike the coordinate formulation, the algebraic Bethe ansatz 
does not require explicit form of the Hamiltonian. The information about the Lax operator 
and the i?-matrix is more important for this method. However using the definitions ( |5.5| ) 
and ([Oj) we can construct Hamiltonian of the model from the Lax operator. For example, 
the XXZ chain with Lax operator fl2,4|) ( with a shift A — > A + -) satisfies the condition 
Aii(0) = Pai = I + cr • a, where P a i is the permutation operator with the property P 2 = 1 
and P a iL' ai+1 (0) = L' ii+l (0)P a i. Using this property we derive from (|3.1[) 

r(0) = tr a (PajPaj+l ■ ■ ■ Paj-i) = (Pj j+ i ■ ■ ■ P 3J ~i)tr a {P aj ) (5.12) 

for any j, applying freedom of cyclic rotation of matrices under the trace. Taking derivative 
with respect to A in ( |3.1[ ) we similarly get 

N N 

r'(0) = tr a J2 (P aj K j+ M ■ ■ ■ Paj-l) = E(4i+i(°) • • • Pn-i)tr a (P aj ), (5.13) 

j=l 3=1 

where we have assumed the periodic boundary condition: L aX+ j = L a j. Defining now from 
H = c C\ = lnr(A)u=o = c r/ (0) r_1 (0)) c = const, we get the Hamiltonian from 
( fH2l) and (|5~ll|) as 

N 

H = cE(^'+i(0)^+i) (5-14) 

3=1 
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with only 2-neighbor interactions, due to cancelation of all other nonlocal factors and the 
property = P. Calculating V from the Lax operator (2.4) and inserting in ( |5.14| ) we 



finally get the explicit form of H as in (2.4) for the XX Z chain. 



6 Algebraic structures in ultralocal and nonultralocal inte- 
grable systems 

We have seen that the algebraic property in integrable systems is dictated by the QYBE. 



However since the QYBE for the nonultralocal systems is different from (3^2), the related 
algebraic structure can not be described by the above procedure. As a result progress 
in this field is not much satisfactory compared to the above described ultralocal theory, 
though many important models like Nonabelian Toda chain, quantum KdV and modified 
KdV model, nonlinear cr-model etc. belong to this class. Nevertheless, it is possible now 



to describe such models to considerable extent by some extensions of the QYBE [^|] [23]. 

However for better understanding of the algebraic structure of the nonultralocal models 
let us take a closer look at that of the ultralocal Lax operators themselves. The underlying 
quadratic quantum algebra Q4.2[ ), as mentioned before, exhibits Hopf algebra property. The 
most prominent characteristic of it is the coproduct structure given by 



A{s 3 ) = s 3 ®/+/®s 3 , A(S + ) = S + ®4q- s3 +c+q s3 ®S + , A(S~) = S^c^ 3 +c^q s %S~ 

(6.1) 

Note that the = 1 case recovers the well known result for the standard quantum algebra. 
This means that if = (g) I and = I ® satisfy the quantum algebra separately, 
then their tensor product A(S ,± ) given by (|6.1|) should also satisfy the same algebra. This 



Hopf algebraic property induces the crucial transition from the local QYBE (3.2) to its 
tensor product given by the global equation ([^4|) as discussed above. Note that thanks to 
the property of the generators of the algebra: [Sf, Sj] = 0, i ^ j the ultralocality of the 
Lax operators is achieved. However for nonultralocal models with the underlying braided 



algebra having a different multiplication rule though a similar Hopf algebra property [21|, 
we get [Lu,L2j] / requiring generalisation of the QYBE. 

Such generalised QYBE for nonultralocal systems with the inclusion of braiding matrices 
Z (nearest neighbour braiding) and Z (nonnearest neighbour braiding ) may be given by 

R\ 2 {u - v)Z2i(u,v)Lij(u)Z2i(u,v)L 2 j(v) = Z^ 2 1 (v.u)L 2 j{v)Zi 2 (v,u)Lij(u)Ri2(u - v). 

(6.2) 

In addition, this must be complemented by the braiding relations 

L 2j+ i(v)Z 2 ~ 1 1 (u,v)L lj (u) = Z 21 1 (u 1 v)Li j (u)Z 21 (u,v)L 2j+1 (v)Z 21 l (u,v) (6.3) 
at nearest neighbour points and 

^k{v)Zn(u,v)Lxj(u) = Z^-(u.v)L lj (u)Z 2 i(u,v)L 2 k(v)Z2i(u,v) (6.4) 



12 



with k > j + 1 answering for the nonnearest neighbours. Note that along with the usual 
quantum R\2(u — w)-matrix like fl3.5| ) additional Z\2, Z\2 matrices appear, which can 
be (in-)dependent of the spectral parameters and satisfy a system of Yang-Baxter type 
relations [B3|. 



7 Generation of nonultralocal models from the braided 
QYBE 

Unlike ultralocal models due to the appearance of Z matrices in the braided QYBE relations 
one faces initial difficulty in trace factorisation. Nevertheless, in nonultralocal cases one can 
mostly bypass this problem by introducing a K(u) matrix and defining t(u) = tr(K(u)T(u)) 
as commuting matrices |3^, ^] for establishing the quantum integrability for nonultralocal 
models. Though a well-framed theory for such systems is yet to be achieved, one can derive 
the basic equations for a series of nonultralocal models in a rather systematic way from the 
general relations ( |6.2| - |6.4D by paricular explicit choices of Z, 2 and i?-matrices (23|, [2£| . The 
models which can be covered through this scheme are 
I. Nonultralocal models with rational R-matrix 
1. Nonabelian Toda chain |24| 

Z = 1, Z = I + ih(e22 <S> ei 2 ) <8> 7r. 



2. Nonultralocal quantum mapping [25| 

2 = 1 and #12(1*2) = 1 + ^ J2a _1 e Na <8> e aN . 

3. Super 'symmetric models 

Z = 2 = Y.fla.pgo.p, where r) a/3 = e aa <g> epp and g = (-l) a/3 with super symmetric 
grading a. 

4. Anyonic type SUSY model 

Z = 2 = J2 r]apg a p, with g a p = e lda ^ ' . 

5. Kundu-Eckhaus equation [36| 

Classically integrable NLS equation with 5th power nonlinearity 

iipi + ^xx + n{ip^)i> + ^(VV) 2 ^ + 2^(VV)^ = 0, (7.1) 

as a quantum model involves anyonic type fields: rp n ip m = e '.ipmipm n > m; [ijj m , tp^ ] = 1. 
The choice Z = 1, Z = diag(e %e , 1, 1, e l8 ) constructs the braided QYBE, The trace factori- 
sation problem has not been solved. 

II. Nonultralocal models with trigonometric R-matrix 



1. Current algebra in WZWN model [19| 



Z = 1 and R12 = Z\2 = R q i2, where R^ is the A — > ±00 limit of the trigonometric 
i?(A)-matrix 



2. Coulomb gas picture of CFT [26| 



2 = 1, Z 12 = q -Y,i H i® H i and ,Ri 2 = R+ 12 . 
3. Integrable model on moduli space \27\ 
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Z = Z\2 = Rq and A-dependent i2(A)-matrix 
4. Quantum mKdV model pgj 

Z = 1, Z12 = ^21 = q~? a ® CT , and the A-dependent trigonometric Z2(A)-matrix. 

Other models of nonultralocal class are the well known Calogero- Sutherland (CS) and 
Haldane-Shastry (HS) models with interesting long-range interactions. Though spin CS 
model exhibits many fascinating features [^] and its integrability formulation through 
braided QYBE for both HS and CS models has not yet been achieved. 



8 Connections with Stat Mech and CFT 



As it was mentioned above the unifying feature of the QIS is present in its relationships with 
various other branches of physics and mathematics, where the knowledge and techniques 
of the QIS often proved to be extremely helpful. Among these various connections those 
with the Stat Mech problems and the CFT models seem to be the most interesting. 

The (1 + 1) dimensional quantum systems are intimately linked with 2 dimensional 
classical statistical models. Moreover the notions of integrability are equivalent in both 
these cases and are through the Yang-Baxter relations. In integrable statistical systems 



both QYBE (^2) and the YBE ( |3.3| ) become the same and leads similarly to the commuting 
transfer matrices r(A). A classical statistical vertex model may be given by N x M lattice 
points connected by the bonds assigned with +ve (-ve) signs or equivalently, with right, up 
(left, down) arrows in a random way (see Fig. 1). 



n 



4 iy. 



v. 



s — - 



Figure 1 Classical 2 -dimensional vertex model is related with (1 + 1) -quantum system. 
The row-to row transfer matrix t depends only on vertical space indices a, (5 and made out 
of multiplying -matrices with elements R mn, which represent the Boltzmann weights 
at vertex point i. 
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Setting the corresponding Boltzmann weights for possible arrangements as matrix el- 

(i) 

ements of a 4 x 4-matrix, we get the R u ; -matrix with crucial dependence on spectral 
parameter A. The YBE ( |3.3| ) restricts the solution of the i?-matrix to integrable models. 
Imposing extra symmetries on the i?-matrix like the charge conserving: i?^z 0> > only 
when k + I = i + j and reversing = R~l~~?i symmetry we get the Boltzmann weights of 
the 6-vertex model constituting the elements of the i?-matrix: 

R$+ = RZZ = a(A), R+Z = RZX = b{\), 

Rtl = Rlt = c(A), 

where f = /, f = /i- Using an overall normalisation one immediately recognises this as our 
familiar i?-matrix (]3.5|), which represents also the Lax operator of the anisotropic quantum 
spin-i chain fl2,4|) . 

The configuration probability for a string of iV-lattice sites in a row may be given 
by the transfer matrix rjy(a, fj) = tr{J\f R^) in a similar way as in QIS. The free 
parameter A may be linked with the temperature T = of the statistical model as 
^ _ const, rpj ie direction along N in statistical system may be considered as the space 
direction of the quantum model, while the M direction is like the time direction in the 
quantum system (see Fig. 1), which again may be replaced in statistical systems by the 
temperature. Therefore, the finite-temperature behaviour can be studied by taking suit- 
ably the limit of M [Oj. Recalling that the Hamiltonian of the quantum model may be 
given by H = C\ = — tn -^tn (X) \x=o, we can establish the important connection 

exp(—f3H) = limjvf^oo (t® td ^ , where t® td = T7v(0) _1 ttv(^) is the diagonal to diago- 
nal transfer matrix. On the other hand the partition function Zn.m(^) = tr (r® TD '(A)^ 
Interchanging N, M amounts to the rotation by 90° and effected by changing the spectral 
parameter A — > r/ — A due to the crossing symmetry of the model. This gives therefore 
Zm.nW = tr \T^[ TD (r] - A) J yielding the free energy / = -± lim M ^oo A M (A), where 
Am (A) is the largest eigenvalue of t^ d {j] — A). Therefore using the eigenvalue expressions 



( |5.S| ) and expanding it in powers of /3 one can find the finite-temperature corrections to 
the thermodynamical quantities just as in the case of expansion in large N for finding the 
finite-size corrections as shown below. 

Similar deep interrelation exists also between the QIST and the CFT models, revealed 
first perhaps by Zamolodchikov [|l2| by showing that, if CFT is perturbed through relevant 
perturbation and the system goes away from criticality, it might generate hierarchies of 
integrable systems. For example c = \ CFT perturbed by the field a = 4>(i,2) as H = 
Hi+ha J a(x)d 2 x, represents in fact the Ising model at T = T c with nonvanishing magnetic 

2 

field h. Similarly the WZWN model perturbed by the operator ^(^3) generates integrable 
restricted sine-Gordon (RSG) model. Under such perturbations the trace of the tress tensor, 
unlike pure CFT becomes nonvanishing and generates infinite series of integrals of motion 
associated with the integrable systems. 

Another practical application of this relationship is to extract important information 
about the underlying CFT in the scaling limit of the integrable lattice models. Interestingly, 
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from the finite size correction of the Bethe ansatz solutions, one can determine the CFT 



characteristics like the central charge and the conformal dimensions [13]. For example, one 
may analyse the finite size effect of the Bethe ansatz solutions of the 6- vertex model ( with 

i — 

a seam given by k). Considering the coupling parameter q = e "+ 1 , one obtains from the 
Bethe solution at the large N limit the expression 

for the ground state energy and 

2-7T ~ 1 2tt ~ 1 

E m - E = — (A + A) + O(^) P m - P = —(A - A) + O(^) 

for the excited states. Here A, A are conformal weights of unitary minimal models and 



c = 1 — ^7^j?r) i v = 2,3,... is the central charge of the corresponding conformal field 
theory. 

9 Concluding remarks 

The variety of models within the quantum integrable systems are linked together by an 
unifying algebraic structure induced by the quantum Yang-Baxter equation. Such struc- 
tures not only systematically generate different classes of integrable models but also gives 
an unifying approach in their Bethe ansatz solutions. The theory of nonultralocal models, 
though not so successful allows to derive the models from few basic relations, which are 
braided extensions of the Yang-Baxter equations. The relationship of QIS with stat mech 
and CFT gives practical results in finding finite-temperature and finite size corrections of 
the quantum model as well as its underlying conformal properties. 



Recently found link of the QIS with diverse subjects like link and knot polynomials [15], 



reaction-diffusion processes p9fl , Seiberg-Witten model [ 40 1 etc. are also becoming more 
and more important. 
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